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ABSTRACT 

We have analyzed the redshift-dependent fraction of galactic bars over 0.2<z<0.84 in 2,157 luminous 
face-on spiral galaxies from the COSMOS 2-square degree field. Our sample is an order of magnitude 
larger than that used in any previous investigation, and is based on substantially deeper imaging data 
than that available from earlier wide-area studies of high-redshift galaxy morphology. We find that the 
fraction of barred spirals declines rapidly with redshift. Whereas in the local Universe about 65% of 
luminous spiral galaxies contain bars (SB+SAB), at z ~ 0.84 this fraction drops to about 20%. Over 
this redshift range the fraction of strong (SB) bars drops from about 30% to under 10%. It is clear 
that when the Universe was half its present age, the census of galaxies on the Hubble sequence was 
fundamentally different from that of the present day. A major clue to understanding this phenomenon 
has also emerged from our analysis, which shows that the bar fraction in spiral galaxies is a strong 
function of stellar mass, integrated color and bulge prominence. The bar fraction in very massive, 
luminous spirals is about constant out to z^0.84 whereas for the low mass, blue spirals it declines 
significantly with redshift beyond z=0.3. There is also a slight preference for bars in bulge dominated 
systems at high redshifts which may be an important clue towards the co-evolution of bars, bulges and 
black holes. Our results thus have important ramifications for the processes responsible for galactic 
downsizing, suggesting that massive galaxies matured early in a dynamical sense, and not just as a 
result of the regulation of their star formation rate. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: spiral — galaxies: struc- 
ture — galaxies: general 



1. INTRODUCTION 

How, when and at what rate did the Hubble sequence 
form? This question is central to the field of galaxy for- 
mation and evolution. We examine it by measuring the 
evolution of the bar fraction with redshift using the 2- 
square degree Cosmic Evolution Survey (COSMOS). In 
nearly all simulations, the formation timescale for a bar 
is rapid once the necessary conditions (a massive, dy- 
namically cold and rotationally-supported disk) are met. 
Therefore the redshift evolution of the bar fraction is a 
fundamental probe of the evolutionary history of disk 
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galaxies. 

The bar fraction is defined simply as: 

number of barred spirals 

number of all spirals 

In the local Universe the value of fbar is quite well es- 
tablished. When only strongly bar red 12 galaxies (SB) are 
counted, the RSA, RC3 and UGC (ISandage fc Tammaim] 
119871: iDe Vaucouleurdll99lUNilsonlll973f) all give values 
of f bar = 0.25 - 0.3. When ovally distorted (SAB) are 
also counted the situation becomes a little less clear-cut, 
because, unlike the RC3, the UGC and RSA do not at- 
tempt to carefully compile an inventory of such galax- 
ies. If ovally distorted systems in the RC3 are included 
in the computation of fbar then the local bar fraction 
rises to fb ar ~ 0.6. This result is in good agreement 
with recent infrared studies whic h have measured the 

0.65 (lEskridgeet aI1 l200C: 
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Whvte et al. 2002; Mcncndcz-Dcl mestre et all 
Marinova fc Jogegl2007h . In the infrared, a majority of 
the S AB galaxies are classifi ed as strongly barred SB sys- 
tems (lEskridge et al.ll2000D. As noted b y lEskridge et all 
(2000) and iMenendez-Delmestre etal] (|2007[ ). the over- 
all bar fraction is the same in the infrared and the optical 
(although there is a small number of cases where bars are 
unveiled at infrared wavelengths) . This is not surprising, 
because bars are primarily stellar structures whose vis- 
ibility only declines sharply at ultraviolet wavelengths, 

12 Bars that are highly elliptical and have rectangular isophotes 
are classified as strongly barred (SB) galaxies whereas those with 
more oval shapes are classified as SAB or ovally distorted galaxies 



2 



Sheth et al. 



short wards of the Balmer break (see also ^A. 1[) . We con- 
clude that the consensus value of the local barred fraction 
is fbar ~ 0.3 for strongly barred systems, and fbar ~ 0.65 
for all barred galaxies, and that these values are so well- 
known that they have not changed significantly in over 
four decades. 

In sharp contrast with the rapid and stable con- 
sensus reached on the local bar fraction, attempts to 
measure the bar fraction at high redshift have proven 
difficult. The earliest analyses of the bar fraction 
in the Hubble Deep Fields (HD Fs) found a dramatic 
paucity of barred spirals at z>0. 5 (Abraham et al.1119961 
Ivan den Bergh et alJll996t lAbraham et al.lll999h . These 
authors concluded that at lookback times greater than 
5 Gyr disks were either dark matter dominated or dy- 
namically too hot (perhaps due to the increased merg- 
ing activity) to host bars. However, the small volume 
probed by the HDFs (only thirty bright, face-on spi- 
ral galaxies between < z < 1) led to concerns that 
the bar frac tion at high r edshif t may not be adequately 
measured. ISheth et all (|2003f ) investigated whether a 
significant n umber o f bars could have been missed, as 
sugge sted bvlBunkeri (H9 99) using the H-band NICMOS 
HDF. [Sheth et all (|2003rT found four bars and two can- 
didate bars out of 95 galaxies at z>0.7. Overall, the 
fraction of barred spirals in the NICMOS HDF remained 
extremely l o w, as in the optical HDF studies. But 
ISheth et al.1 {2003) noted that their study was limited 
by the coarse NICMOS resolution (0'.'15) such that only 
the largest (and rarest) bars could be identified (bars 
with semi-major axis >5 kpc). When the fraction of 
these large bars at z>0.7 was compared to local samples, 
there was no compelling evidence for a decline in barred 
spirals, but likewise the NICMOS data did not unveil 
any new b ars at low redshifts ; all except one of the four 
bars in the ISheth et alJ (|2003f ) stud y are at z>0.9, where 
k-correction effects are important ( §A.ljl . 

A major advance in spatial resolution was possible 
with the Advanced Camera for Surveys (ACS) whose 
0'.'05 pixels are able to resolve all but the smallest (nu- 
clear, <2 kpc diameter) bars at all redshifts. Using 
ACS data, two studies ([Elmegreen. Elmegreen fc Hirstj 
120041 Uogee et al.1 l2004h found that contrary to the 
previous HDF results, the bar fraction is constant at 
30% over the last 8 Gyr (since z = 1.2). The sam- 
ple siz es, however remained modest in th ese studies 
(186 in [ Elmegree n. Elmegreen fc HirstJ 12004. and 258 in 
Uogee et alJl2004D . 

In this paper we examine in detail the redshift evolu- 
tion of the bar fraction using the unparalleled wide and 
deep 2-square degree COSMOS data set. The plan for 
the paper is as follows: in Section 2 we describe our sam- 
ple selection procedure. The classification methodology 
we have adopted is described in Section 3. Our main re- 
sults are presented in Section 4, before being discussed 
in Section 5. Our conclusions are summarized in Section 
6. An Appendix to this paper provides a detailed anal- 
ysis of possible selection effects at high redshift and a 
discussion of our local calibration sample of 139 galaxies 
fr om the Sloan Digital Sky Survey (SDSS) Data Release 
4 (jAdelman-McCarthv et al.ll2006h . Throughout this pa- 
per we adopt a flat A-dominated cosmology with Hq=70 
km s" 1 Mpc" 1 , VL M = 0.3, and Q A = 0.7. 



2. SAMPLE SELECTION 

An overview of th e COSMOS program is given in 
IScoville et all (|2007a[ ) and details of the H ST observa- 
tions are described in IScoville et~aTl (|2007bf L The COS- 
MOS observations are undertaken in the F814W fil- 
ter ('I-band') and reach a depth of Iab > 27 (10cr). 
The photometric catalog an d redshift measur e ments used 
in this pape r are given in Mobash er et al.l (|2007h and 
iCapak et al.1 (|200l . 

The most important and difficult step in studying the 
evolution of galactic structures is choosing comparable 
samples at different lookback times. For nearby galaxies, 
multi- waveband data with sufficient spatial resolution are 
available. Therefore galaxy properties (e.g., Hubble type, 
inclination, distance) are accurately known for nearby 
spirals. The underlying disk is also well-imaged and mul- 
tiple techniques for identifying a bar may be employed. 
However, for high redshift galaxies, the situation is more 
complicated. Here we summarize the steps we have taken 
to overcome these problems. 

We choose all galaxies brighter than L v with an em- 
piric ally dete r mined luminosity evolution of 1 magnitude 
from ICapakl f2003). such that M^> = -21.7 at z=0.9 
(|Capakll2003f ). This criterion is specifically targeted for 
choosing galaxies from the same portion of the galaxy 
luminosity distribution at all redshifts. As we shall see 
later, a no-evolution luminosity model would have only 
steepened the overall decline in the bar fraction, fur- 
ther strengthening the results presented in this paper. 
For this luminosity criterion, at a redshift of z=0.9, the 
sample is complete for all galaxies with a half-light ra- 
dius s maller than 10 kpc (see Figure 6, IScoville et al.1 
l2007bf ). The number of galaxies with a half light radius 
larger than 10 kpc is extremely lo w (e.g., Figure 10 in 
ISargent et ai1l2007l or Figure 10 in iBarden etalll2005h 
and thus our sample is essentially complete. 

We eliminate all elliptical and lenticular galaxies based 
on the galaxy's spectral energy distribution (SED) type 
Tphot, choosing all galaxies with T phot > 2. T phot is the 
best fit spec tral template ordered by the 4000A break 
strength (see lMobasher et aT1l2007t ). The Tphot sequence 
can be thought of as a photometric Hubble type go- 
ing from the reddest early type (T p ^ ot =l) galaxies to 
bluest late type (T p ^ ot =6) galaxi es. Types 1 through 4 
are defined by the templates from lColeman et alJ (1980) 
and correspond to Elliptical, Sbc, Scd, and Irregular 
H ubble types re s pectiv ely. Types 5 and 6 correspond 
to iKinnev et al.l (|1996l ) type SB3 and SB2 respectively 
which are local star-burst galaxies with little or no ex- 
tinction. Typical uncertainties in Tphot are ±0.2. 

The Tphot values are robust descriptors of galaxies to 
z^l.2 as confirmed by a comparison of the photometric 
and spectroscopic redshift s for over eight hundred galax- 
ies (|Mobasher et al.ll2007l ). Comparison of galaxy types 
based on morphological parameters such as Gini and 
asymmetry and the T v ho t type shows that T v hot > 2 se- 
lects all spirals to z— 1.2 (|Capak et al.ll2007t lllbert et all 
12001 . With increasing redshift there is an increase in the 
population of blue ellipticals and blue-bulge dominated 
spirals; however these galaxies are only signifi cant at the 
faint end of the luminosity function (e.g ., lllbert etldl 
l2006UCapak et al.ll2007K lllbert et al1l2007f ). For instance, 
at Mb <-20, the blue-bulge dominated spiral population 
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is les s than 1% of the total disk population (jllbert et al.l 

2006) . At z~0.8, amongst elliptical galaxies (identified 
by the Gini and asymmetry parameters), the contribu- 
tion of blue ellipticals to the total volume density is sig- 
nificant at the very low mass end (~10 9 Mq). For the 
high luminosity (massive) galaxies, which are studied in 
this sample, the contribution from blue ellipticals is no 
more than a few percent. We do find a small fraction of 
compact systems with blue colors which we identify and 
reject from our analysis. 

We impose a redshift cutoff of z = 0.835 so that the 
F814W filter does not probe bluer than rest-frame g- 
band. The rationale for probing no bluer than rest-frame 
g-band is described in detail in ^A.ll of the Appendix, 
but we also note here that at this redshift the angu- 
lar diameter of a 0.05" ACS pixel subtends a physical 
scale of 0.4 kpc. At this resolution we expect to de- 
tect easily the entire po pulation of bars in nearby spiral 
galaxies (see Figure 3 of Shcth ct al. 2003, or Figure 7 of 
IMenendez-Delmestre et aT1l2007t ). We note that all our 
galaxies subtend at least 10 pixels. The smallest pet- 
rosian radius in our sample is 5.2 pixels. 

Imposing these cuts reduces our sample to 3886 spirals. 
However, galactic structures like bars are difficult to iden- 
tify and quantify in edge-on galaxies. Therefore, wc fur- 
ther eliminate all galaxies with inclinations i > 65°, the 
same limiting inclination use d in studies of the bar frac- 
tion in nearby spirals (e.g., IMenendez-Delmestre et al.l 

2007) . Inclination values for our sample are determined 
from the axial ratio of the galactic disk, which is identi- 
fied using a two dimensional decomposition with a bulge 
and e xponential compo nent for each galaxy using GAL- 
FIT (Pen g et al.ll2002l ). The inclination cut eliminates 
an additional 986 objects, leaving 2900 galaxies. 

Finally, we discard all galaxies that were obviously 
merging, or too irregular or peculiar to be fit with el- 
lipses. This eliminates 743 galaxies, leaving a final sam- 
ple of 2,157 spiral galaxies which were then classified as 
barred or unbarred. The importance of eliminating pe- 
culiar objects from our analysis is investigated further in 
§A.2I of the Appendix, but we note here that including 
these objects would not have changed any of our conclu- 
sions. 

All galaxies are detected in at least thirteen, sometimes 
all sixteen available photometric bands. The galaxy lu- 
minosities, masses and c olors are measured fro m the pho- 
tometric redshift code (Mobashcr ct al. 2007) which si- 
multaneously solves for observed galaxy flux, redshift, 
galaxy type, and extinction intrinsic to the galaxy. The 
galaxy type combined with the flux-normalization yields 
rest frame, extinction corrected, luminosities and colors 
for each object. The rest frame color and luminosity 
is the n converted into stellar mass using the Bel l et all 
(2005) relation. Errors in the mass estimate are within 
a factor of 3 due to systematics. With our deep multi- 
wavelength photometry these parameters are robust and 
non-degenerate at z<1.2. Over this redshift range we 
have at least two data points red-ward of the 4000 A 
break (z + and K s at z=1.2) and two points blue- ward 
(u* and Bj at z=0.3), which allows us to break the 
4000A break strength-extinction degeneracy. More de- 
tailed SED fits using spectroscopic redshifts yield sim- 
ilar results for the galaxy luminosity, mass and color 



(|Mobasher et al.ll2007l : lllbert et al.ll2007t ). 

3. BAR CLASSIFICATION METHODOLOGY 

We identified bars using two methods and the results 
were cross-checked for consistency. We use the same 
methods for the local SDSS calibration sample (see lA.lj) 
and the COSMOS data to reduce biases that can be in- 
duced by the use of different classification methods. 

Our first method was to use the ellipse fitting 
technique that has been used widely by multiple 
studies of both nearby and hi gh r ed shift galax- 
ies (e. g.. iRegan fc ElmegreerJ 119971: ISheth et all 



20051 iKnapen .et al .1 |2000 | : ISheth .et alJ 120021: 



Laine et al.l 120021: ISheth et alJ [20031 : Uoeee et~aTl 12004 



Menendez-Dehnestre et al.l 120071 : iMarinova fc Joged 
2007IL For a detailed discussion of the ellipse fitting 



m ethod and the classification s cheme we refer the reader 
to IMenendez-Delmestre et ail ((2007). Briefly, bars are 
identified from a dual signature of the ellipticity and 
position angle profiles. We require that the ellipticity 
increases monotonically, exceeds 0.2, and then drops 
by at least 5e >0.1. The position angle profile should 
show a relatively constant position angle over the bar 
region and change by at least 10° after the bar as the 
isophotes enter the disk. In general the fitting procedure 
is extremely robust and fits ellipses to the la noise level 
in the images - in the case of COSMOS this means we 
are able to fit galaxies to the outermost disk isophotes 
(at fi=24.5 mag arcsec -2 , see Figure 9) in the highest 
redshift bins. In a small number of cases, the ellipse 
fitting method can miss an exist ing bar (see Figure 3 in 
IMenendez-Delmestre et all [2007). This usually happens 
when the position angle of the disk and bar are aligned 
or when the ellipticity drop is moderated by open spiral 
arms. 

In addition to the ellipticity and position angle pro- 
files, we classified the bars further into strong bars and 
intermediate bars by visually examining the isophote 
shapes. The strong bars are those with an ellipticity 
greater than 0.4 to be consistent with previous work 
in this field (e.g., Uogee et"ai1l2004l L However we note 
that the definition is ar bitrary because there is a con- 
tinuum of bar strengths dBlock et al.l 2002; Whvte et al.l 
120021 iButa. Laurikainen fc Said 120041 iButa et al.l 120051 : 
IMenendez-Delmestre et aill2007t ). We refer to this first 
method as the 'ElPa' classification method for the re- 
mainder of the paper. 

Our second method for classifying galaxies was vi- 
sual identification of galaxies into categories of SA 
(unbarred), SAB, and SB by one of the authors 
(DME). Not all galaxies were classified into these 
three neat categories. DME also classified galax- 
ies as edge-on, clump-cluster galaxy, or compact 
(spher oidal) galaxy following her work o n the Ultra Deep 
Field ([Elmegreen fc Elmegreenl l2005at l Elmegreen et al. | 
2005bjjcj; lElmegreen fc Elmegreenll2006l : lElmegreen et al.l 
2007). Galaxies that were visually classified as edge-on, 
clump-cluster, or compact are not included in the visual 
classification results. The total number of galaxies clas- 
sified visually into barred and unbarred spirals is 1,705. 

The classification by DME was cross-checked for five 
hundred galaxies by another author (KS), and also cross- 
checked against the ellipse fitting method. The cross 
checks between the two methods finds agreement for 85% 
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of the sample. In general we found that we classified more 
galaxies as barred by eye than using the ellipse fitting- 
profile. This is as expected because there are particular 
morphologies where the ellipticity and position angle sig- 
nature can be masked by the relative orientation of the 
bar and disk, and the pitch angle of the spiral arms. As 
noted ea rlier, a detailed discussion of such cases can be 
found in iMenendez-Delmestre et all (|2007l ). Only in 5% 
of the cases is there a gross mismatch where one method 
differed from the other by two classes, i.e. bar in one 
and spiral in the other. This generally occurred for very 
faint or small nuclei where the ellipse fitting method has 
problems fitting the isophotes. In the remaining 10% of 
the cases, the methods agreed to within one class, i.e. 
intermediate bar in one, and a spiral or strong bar in the 
other. The cross checks were internally consistent at all 
redshifts. 

4. RESULTS 

4.1. The Declining Bar Fraction 

Figure Q] shows the evolution of the bar fraction in 
COSMOS as a function of redshift in six equal redshift 
bins (<5z = 0.117), starting at z=0.14 at a lookback time 
of 1.8 Gyr. The redshift bins correspond to lookback 
times of 1.8-3.0, 3.0-4.1, 4.1-5.0, 5.0-5.8, 5.8-6.5, and 
6.5-7.1 Gyr respectively. The two rows show the bar 
fractions measured from the two classification methods 
described above The left column shows the total 

bar fractions (strong bars + oval bars) and right col- 
umn shows the strong-bar fraction (SB). For each data 
point we list the number of bars and the total number of 
galaxies classified in each bin. The error bars reflect the 
statistical uncertainty in the fraction and are calculated 

1/2 

from the expression (/ [1 — /] /N) 1 for fraction / and 
number of galaxies N . 

We find that the bar fraction (for all galaxy luminosi- 
ties combined) has evolved dramatically over the last 
7 Gyr. At z=0.84, the total bar fraction using the 
ElPa classification method is f bar = 0.22±0.02 (0.31±0.02 
for visual classification method), one-third (one-half) its 
present-day value. We see the same trend when we 
consider only the strong bar fraction: fg^ evolves from 
0.27±0.05 (0.35±0.05) at z=0.0 to 0.09±0.01 (0.17±0.02) 
at z=0.835. 

When combined with the z=0.0 data point from our 
analysis of a local SDSS sample (Table 1, Figure O & 
£|A.1[) . we also find that the evolutionary trend is weaker 
in the first three bins at z < 0.4. Within the error bars, 
the data at these redshifts are consistent with a roughly 
constant bar fraction (f Dor =0.6, fss = 0.3). These results 
are summarized in Table [TJ 

4.2. Bar Fraction as a Function of Galaxy Mass & 
Luminosity 

Figures [2] and [3] show %ar versus the absolute lumi- 
nosity and mass of the disk, respectively, in the three 
redshift bins from z=0.14 to z=0.84. We find that in the 
highest redshift bin, galaxies with masses log M (M Q )> 
10.9 and luminosities My < -23.5 have f oar ~0.5, which 
is about the local value. In contrast, the low mass (log M 
< 10.5) and low luminosity (My > -22.5) galaxies have 
{bar < 0.2 at high redshift. The same trend is seen for 
strong bars, tsB- At low redshifts, the bar fraction is 
roughly equal for all luminosities. 



This trend is not due to incompleteness in the sam- 
ple. We establish the completeness of our sample by 
measuring the mass limit based on our selection criteria. 
Since we choose galaxies based on a luminosity cutoff 
and galaxy colors, the mass completeness is most likely 
to be an issue for the reddest systems at the highest 
redshift. For our luminosity cutoff and Tphot criteria, 
our sample is complete for galaxies with masses greater 
than 3-4xl0 10 Af Q at z=0.9 for the reddest (T phot =2, 

rest-frame Am 9 r > 0.56) galaxies. Obviously, for the 

bluest systems (e.g., T p hot = 6), our sample is complete 
to 0.9-1 xlO M Mr. The se values a r e calc ulated from 
the iMarastonl (|2005f ) and iBell et alJ (|2005| ) models re- 
spectively. Note that at z=0.6, the mass limit for com- 
pleteness in the sample is lowered by another ^25%. Our 
lowest mass data point in the 0.6<z<0.84 bin (the high- 
est redshift bin) in Figure [3] is for galaxies with masses 
between 3-4x 10 10 Mq. It is therefore free from the possi- 
ble mass selection bias. The data points at lower masses 
and lower redshifts are also computed from a complete 
sample of masses for a given bin. Thus we conclude that 
the observed strong correlation between the bar fraction 
and mass in the highest redshift bin is a robust result. 

The most important result in these figures is that in 
the highest redshift bins in this study, a majority of the 
most massive and luminous systems are barred. There is 
little evolution in the bar fraction with redshift in these 
systems. Since bars form in massive, dynamically cold 
and rotationally supported galaxies, the high bar fraction 
indicates that the most massive systems are already "ma- 
ture" enough to host bars. This agrees with the analysis 
of the evolutio n of the size function of disk galaxies of 
several studies (Sargent et al. 2007; Ravindranat h et "all 
12004 iBarden et al.ll2005HSheth et al.ll2007bh . which find 
that large disks are already in place by z=l and little or 
no evolution in disk sizes from z~l to the present epoch. 
Conversely the low bar fraction in the lower luminosity, 
lower mass systems indicates that these systems are ei- 
ther dynamically hot, not rotationally supported and/or 
have not accreted sufficient mass to host bars. Merging 
activity, which is more common at higher redshifts, is 
also likely to affect the less massive systems more severely 
and may be responsible for heating them up more than 
high mass systems. Bar formation may be delayed in 
these hot disks if they are embedded in a massive dark 
matter halo. Although the exact nature of these disks is 
not yet well-known, there is an indic ation that later typ e 
systems may be dynamically hotter (jKassin et al.l [2007). 
We consider these points further in following sections. 

4.3. Bar Fraction as a Function of Galaxy Color & 
Bulge Luminosity 

Figure 2] shows how % ar varies with galaxy SED type 
(Tphot) and redshift. At low redshift, f oar is indepen- 
dent of Tphot, and at high redshift, fbar decreases from 
early (Tp^t <3) to late types (T phot >3). Similarly, f bar 
decreases with redshift more strongly for the late types 
than the early types. This latter trend is consistent with 
the previous result that the bar fraction changes with 
redshift primarily for the low mass galaxies, which tend 
to have late SED types. 

Finally we consider how the bar fraction varies as a 
function of the bulge light in galaxies. Figure [5] shows 
fbar versus the fraction of bulge luminosity in a galaxy 
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Fig. 1. — Evolution of the bar fraction as a function of redshift in equal bins from z=0.0 to z=0.84, out to a look back time of 7 Gyr. 
The bar fraction drops from 65% in the local Universe to about 20% at z~0.84. The fraction of strong bars (SB) drops from about 30% 
to under 10%. The top row shows the results from the visual classification and the bottom row shows the results based on classification 
using the ellipticity and position angle profiles. The left column shows the bar fraction for all galaxies classified as bars, whereas the right 
column shows the same only for the strong bars. The error bars are calculated as y/(f(l — f)/Ni,i n , where f is the fraction of galaxies, and 
N is the number of galaxies in a given category. The numbers above each data point show the total number of bars (or strong bars) / total 
number of galaxies in the bin. 



for different redshift bins. The bulge magnitude is calcu- 
lated from fitting each g alaxy with a Sers ic + exponential 
profile using GALFIT (|Peng et al.ll2002fh The x-axis in 
this figure is the difference between the bulge magnitude 
measured from the GALFIT fitting and the total (disk 
+ bulge) apparent magnitude. We note that the relative 
calibration across redshift bins should be treated with 
caution because we are not correcting for k-correction ef- 
fects that are known to affect two dimensional decompo- 
sition of galaxies. Within a given redshift bin, however, 
the bulge contribution measurement should be robust ex- 
cept for one important caveat. The fitting algorithm is 
not designed to decompose a bar separately. As a result 
the bar light is likely to be split between the exponen- 



tial and Sersic components. If the light profile of a bar 
is exponential, as it is in later Hubble type galaxies lo- 
cally, the majority of that light is likely to be part of the 
exponential component. On the other hand if the bar is 
relatively short and not highly elliptical, its light is likely 
to be added to the Sersic component. The detailed de- 
composition of the bulge+bar+disk will require a more 
sophisticated approach, which is beyond the scope of this 
paper. 

Keeping the above caveat in mind, we find that ff, ar is 
slightly higher for galaxies that are "bulge dominated" 
compared to galaxies which are not bulge dominated in 
the highest redshift bin. In the lowest redshift bin the 
slight trend disappears and the bar fraction is roughly 



Sheth et al. 



1.0 



Bar Fraction versus Absolute Magnitude 



= 0.6 



0.4 



0.2 



0.4 



0.2 



0.0 



7/11 



-| 1 1 1 1 1 i- 

z=Q.14--0.37 



z=0. 37 — 0. 60 
z=0.60--0.84 



>0/36 



28/43 



27/39 



49/81 



'v^ 44/87 



7/41 



73/126 56/98 

77/181 

--ju^ 91/245 



66/216 " ~f ~ 



- J33/368 

"" % - - _ 82/459 



H — I — | — I — I — I — I — | — I — I — I — I — | — I — I — I — I — | — I — I — I — I — | — I — I — I — I — | — I — I — I — I — | — H 



z = 0.14--0.37 
z = 0.37--0.60 

z=0.60--0.84 



22 



9 22 



34 



„ <► 



24 



44 



32 \ 26 



36 



31 



-24.0 



-23.5 



-23.0 -22.5 -22.0 
Absolute Mag (M_V) 



-21.5 



-21.0 



Fig. 2. — Total bar fraction (top panel) and strong bar fraction (bottom panel) versus the absolute magnitude (My) in three redshift 
bins. There is a strong correlation between the bar fraction and the galaxy luminosity in the highest redshift bin. In that bin, the most 
luminous galaxies (My < -23.5) have f;, ar ~ 0.5 whereas the lowest luminosity (My > -22.5) galaxies have a f;, ar ~ 0.2. With redshift we 
see a strong evolution in the low luminosity sample as they evolve to % ar ~0.6 in the lowest redshift bin. A similar trend is seen in the 
bottom panels for the strong bar fraction. The data points are at the mid-points of bins of <5My=0.5, from -21.0 to -23.5 (a data point is 
skipped if no galaxies are found in a bin) and the errors bars are calculated as before. The first bin is from My = -23.5 to -24.75. The 
data points are slightly offset along the x-axis for each redshift bin and only the number of strong bars is labelled on the bottom panel for 
clarity. This figure should be viewed together with Figure [3] 



constant for all types of bulges, although there are only 
a few galaxies that are not bulge dominated. This cor- 
relation together with the previous correlations (Figures 
[H [31 SI) suggests that the galaxies that are red, luminous 
and massive are also bulge dominated, and in these galax- 
ies the bar fraction does not vary strongly with redshift. 
We discuss the implications of the bar-bulge correlation 



5. DISCUSSION 



5.1. 



Comparison with Previous Studies 

In Figure[6]we plot all of the data from previous studies 
of the evolution of the bar fraction for comparison to the 
COSMOS results. In all cases there is a general decline 
in the bar fraction. However the interpretation of the 
data has been very different amongst these studies as we 
discuss below. 



Our basic result of a decline in the bar fraction is con- 
sistent with the earliest HDF s t udies (jAbraham et al.l 
119961 Ivan den Bergh et al.l [l99l lAbraham et all 119991) 
shown with the red data points in Figure [51 These stud- 
ies reported a striking decline in the fraction of bars at 
z>0.5. But re-analysis of their Figure 4 shows that, in 
fact, at z>0.5, there are 10 barred spirals out of 29 galax- 
ies, and at z>0.8, there are 3 barred spirals out of 11 
galaxies, consistent with the COSMOS results presented 
here. Similarly our result is consistent with the very low 
bar fraction (5-10 %) measured from the NICMOS HDF 
(|Sheth et al.ll2003f) at z>0.7. But it is difficult to com- 
pare this bar fraction to the ours because the NICMOS 
data can only probe the largest bars. Of course the vol- 
ume probed by the HDF studies was too small to allow 
evolution in fj, ar to be probed with much confidence. 

On the other hand, our central result is fundamen- 



COSMOS Bars 



Bar Fraction versus Galaxy Mass 



0.4 



0.2 



0.4 



0.0 



2/5 



29/45 



54/84 'I 7 ' 8 27/43 



\ 51/104,' 

r 

45/1 J i t 



63/12 5 



6/ 



63/169, 



38/84 



25/39 



48/100 
M 



31/204 J 



58/327_ , - J ' 



60/197, 
58/255 , ' i" 



H 1 1 1 1 1 1- 



19 



10 25 



29 



22 " 



13 3r 



H 1- 



12 



10.5 

Log (Moss) 



Fig. 3. — Total bar fraction (top panel) and strong bar fraction (bottom panel) versus total galaxy stellar mass in three redshift bins. 
As expected from Figure [2] there is a strong correlation between the bar fraction and the mass for the highest redshift bin. In that bin, 
galaxies with log M> 10.9 already have f(, ar ~ 0.5, whereas galaxies with log M < 10.5 have fi, ar < 0.2. The bar fraction for the entire 
population evolves with time with the largest change in the lowest mass bin. The same trend is seen in the bottom panel. The lack of 
high mass galaxies in the lower redshift bins is because, even with 2-square degrees, the volume of space observed by COSMOS is small. 
Luminosity and color selection criteria put a limit on the minimum detectable mass - our sample is complete for the points shown here. 
Each point is at the left edge of bins of 51og M=0.15, starting from 10.0 (point skipped if no data is found in a bin). Errors bars are 
calculated as before and last bin is from log M = 10.9 to 11.5. The data points are slightly offset along the x-axis for each redshift bin and 
only the number of strong bars is labelled on the bottom panel for clarity. The uncertainty in the mass measurement is a factor of three. 
This figure should be viewed together with Figure [2] 

galaxies (My <-19.3), which should have resulted in a 
lower bar fraction, their i$B is 50% higher than the COS- 
MOS results. Some possible reasons for the discrepancy 
in the value of fsB may be the different use of incli- 
nation cuts (they used i >60° compared to our cut at 
i > 60°) and selection criteria. We chose our sample 
based on a galaxy luminosity with an evolving luminos- 
ity funct ion, spectral type / color and visual classification, 
whereas [Jogee et all (|2004l ) chose their sample based on 
a fixed (lower) luminosity, range of U-V colors (which is 
similar to our cut in T p hot), and a Sersic parameter from 
a single fit to the galaxies or a concentration index. Given 
the various uncert ainties at hand, w e conclude that the 
data presented by I Jogee et"aT1 ((2004) can be interpreted 
as being consistent with the more significant decline seen 
in our sample. 



tall y different from that reported by I Jogee et alj (|2004D 
and lElmegreen, Elmegreen fc Hirstl (|2004| ); these stud- 
ies have reporte d a constant bar fr action with redshift 
to z~l. In the I Jogee etlrtl (|2004[ ) study, the authors 
classified 258 galaxies as either strongly barred or un- 
barred. Their modest sample size however prevented 
them from studying the fi, ar evolution in detail. Their 
first bin, for example, encompasses our central four red- 
shift bins. Nevertheless we can compare their measure- 
ment of {sb to ours. Over the same redshift range 
in the COSMOS data, we measure a strong bar frac- 
tion, isB = 0.23±0.01 for the visual classification, and 
0.17±0. 01 for the ElPa cla ssification methods. In com- 
parison, [jogeii^2_!^SL re I 22lts a f SB ~ 0.3±0.03 13 . 
Even though the I Jogee et al.1 (|2004h study probes fainter 



13 Errors are not reported in the Jogee et al. (2004) study. We 
measure an error for their data using the bar fraction and total 
number of galaxies reported, in the same way as we measured for 
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Fig. 4. — Total bar fraction (top panel) and strong bar fraction (bottom panel) versus spectral type of a galaxy. In the highest redshift 
bin, the bar fraction decreases monotonically toward later spectral types. In contrast, in the lowest redshift bin the bar fraction is constant 
within the error bars across the Tphot sequence. The data thus show that the majority of the evolution in f(, ar is in the bluer, later SED 
type systems. This is consistent with the results shown in previous figures for high and low mass galaxies, considering the usual correlation 
between mass and Tphot- As before, the data points are slightly offset along the x-axis for each redshift bin and only the number of strong 
bars is labelled on the bottom panel for clarity. 

Elmcgreen, Elmc green k. Hirst] (|2004f ) also reported a 
constant bar fraction to z ~ 1.1 based on an analysis of 
186 background galaxies larger than 10 pixels in diame- 
ter in the multi-color ACS image of the Tadpole galaxy. 
These data points are shown with the the purple triangles 
in Figure [5J The data show a declining bar fraction from 
~ 30% to ~ 15% out to z = 0.8 with 3 — a uncertainty, 
and a rise in the bar fraction from z = 0.8 — 1.1, which is 
beyond the redshift investigated here 14 . Their conclusion 
that the bar fraction is flat on average followed primarily 
from the second rise at z ~ 1; otherwise their fractions 
agree with ours to within statistical uncertainties. Our 
results are also in line with a recent analysis of the Hub- 
ble Ultra Deep Field where t he bar fraction, shown w ith 
the blue triangle, is ~10% ( El megreen et al.l [2005bT ) at 

14 Four galaxies in their Fig. 10 at z ~ 1.8 are incorrect because 
of a photometric redshift error - the corrected data point is shown 
here 



z~l, consistent with the previous HDF studies and the 
values obtained in this paper. 

Although we have attempted to put all the data from 
various studies into context, we emphasize that it is not 
straightforward to make direct comparisons because of 
different selection criteria and bar identification meth- 
ods between these studies. These may be responsible 
for some of the observed differences. The main point to 
note is that in nearly every study, the data have shown a 
decline in the bar fraction, although the interpretations 
of the data have ranged from a constant bar fraction to 
a dramatic paucity of bars at z~l. It is only with the 
COSMOS data set that we are able to robustly quantify 
the decline in the bar fraction and show that the evolu- 
tion is a strong function of the galaxy luminosity, mass, 
color and bulge dominance. 
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Fig. 5. — Bar fraction as a function of the bulge contribution in galaxies in three redshift bins, as before. We find that in the two high 
redshift bins there is a slight preference for bars to be in bulge-dominated systems. The difference is less pronounced in the lowest redshift 
bin. The x-axis is calculated by subtracting the apparent bulge magnitude from the total galaxy magnitude. The bulge magnitude is 
calculated by fitting each galaxy with an exponential and bulge profile using GALFIT as discussed in the text. The x-axis thus indicates 
the fractional contribution of the bulge to the total luminosity of the galaxy with bulge dominated systems to the left on the x-axis. Since 
k-corrections are important in two dimensional decomposition of galaxies, the results across redshift bins are not as robust as those within 
a given redshift bin. The data points are at the midpoints of each T p hot tyP e bin, and as before, only the number of strong bars is listed 
in the bottom panel. 

5.2. Formation of the Hubble Sequence: Assembling the 
Spiral Galaxies 

The declining bar fraction reported in this paper shows 
that at a lookback time of 7 Gyr (z = 0.835) only about 
one-fifth of L* spiral galaxies were barred, which is about 
one-third the present day value. During the following 3 
Gyr (from z ~ 0.8 to z ~ 0.3) the bar fraction increased 
to roughly its present value. Only small changes occurred 
in the last 4 Gyr (z<0.3). 

This evolution can be understood within the frame- 
work of classical bar formation theory. ./V-body simula- 
tions have long suggested that bars form spontaneously 
in galactic disks, usually on relatively fast dynamical 
timescales. There are, however, two ways of slowing this 
down. One is to increase the halo mass fraction within 
the disk radius, and the other is to heat up the disk 



Athanassoula & Scllwood (1986). Self-consistent three- 
dimensional (3D) simulations essentially agree with this, 
although the role of the halo is now understood to be 
more complex, so that the final bar ca n be consider- 
ably s t ronger in cases where it grows slower [Athanassoulal 
(2002, 2003a). Thus, the time it takes for an unbarred 
disk galaxy to become barred can vary widely. In cold, 
disk-dominated cases, the bar forms within a Gyr or less, 
but sufficiently hot disks embedded in very massive ha- 
los can stay unbarred several Gyrs. Such a delay might 
well explain the time evolution of the barred galaxy frac- 
tion shown in Fig. 1. Furthermore, observations show 
that the halo-to-disk mass ratio is higher in low mass, 
low luminosity galaxies than in br ight, massive galax- 
ies (|Bosmall2004l : iKranz et al.ll2003l ) so that bars are ex- 
pected to grow later in the former, as we indeed find 
here. Hence if galactic disks are formed with a variety 
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tional evidence for bar dissolution, ho w ever, is at 
prese nt rather sparse (|Das et al.l l2003t iBlock et"aTI 
2002), while the amount of mass necessary seems 
to be larger than what is currently observ e d for 
central mass concentrations dShen fc Sellwoodl l2004t 



Fig. 6. — Comparison of measurements of our bar fraction with previous studies. f(, ar and fgB measured using the ElPa method for 
the COSMOS data are shown wi th the black fi lled squares and diamonds respectively. The red data circles (13/30 bars in low z bin, 
4/14 in the high-z b in) are from [A braham et al.l j|1999f): t hese do not distinguish between weak and strong bars. Purple triangles are 
bars and twists from Elmcgrccn, Elmcgrccn & Hirst (2004) (10/34, 15/90, 7/31 bars respectively); we summed a djacent bins from their 
data here and corrected the last bin for the incorrect redshifts for four of their galaxies. The green points from Jogee et al. (2003) are 
the strong bar fractions from GEMS; the three points in the two redshift bins are not independent of each other - they are measured for 
~110-175 galaxies, chosen in different ways from the same sample. The horizontal bars show the redshift range over which these data are 
av eraged. Also shown are data from our analysis of a SDSS control sample (diamond - f(, ar , square - fss), and from the 2MASS survey 
by Mcnendcz-Dclmcstrc et al. (2007) (blue diamond is when both ellipticity and position angle signatures are present, and the square also 
includes candidate bars. The vertical dotted line is the limiting redshift for our survey. Within the error bars all the data seem to be in 
agreement (although the Jogee et al. (2004) data point for the lower redshift bin is significantly above the other measurements of fss). 
Contrary to earlier interpretations, it seems that all studies show a general decline in the bar fraction with redshift. It is only with the 
COSMOS data that we are able to analyze this decline in detail. 

of velocity dispersions and a variety of halo-to-disk mass 
ratios, there should be a continuous increase with time in 
the barred fraction, as observed here. The slope of this 
evolution will depend on the distribution of the initial 
disk and halo parameters. On the contrary, if all galaxy 
disks were, in the relevant time-period, similar, then the 
fraction of disk galaxies that are barred would be more 
or less constant with time, or show only a very small 
increase. 

The preceding paragraphs assumed the existence of 
the appropriate set of physical conditions needed to form 
bars from the classical disk instability. Alternatively, 
bar formation may coincide with inner disk growth. 
This, however, would imply considerable growth of the 
inner disk even for z~0.84, which seems inconsistent 
with recent results that s how no evolution in the 
disk scale-len gths at z<l dRavindranath et al.l l2004t 
Barden et all 120051 : ISargent et alJ l2007t ISheth et all 
2007H) . A further point to take into consideration is 
that bars may dissolve when a gaseous component is in- 
cluded in the angular momentum exchange cycle, and/or 
in the presence of a sufficiently massive, centrally con- 
centr at ed object, black hole (iFriedli fe Bend 
1993t iBerentzen et all 119981: iFukuda et all 120001: 
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lAthanassoula, Lambert &: Dehnenl 20051 ). Neverthe- 
less, it is by no means clear that this mechanism is 
unimportant. 

We have so far discussed only isolated galaxies. Let 
us now turn to the effect of interactions and mergings. 
The number of interactions are known to increase dra- 
matically with redshift (e.g.. iKartaltepe et al.ll2007l and 
references therein). Interactions and merging activity are 
most likely to influence (heat up) the less massive galax- 
ies. It is precisely in such galaxies that we see signif- 
icantly lower bar fractions compared to the high mass 
galaxies at the highest redshifts (Figure [3l [4]) . Although 
indirect, there is observational evidence that later type 
and less massive systems are dyna mically hotte r. The 
top row of panels in Figure 1 of iKassin et al.l (|2007f ) 
clearly shows that late-type spirals and irregulars have 
larger disordered motions compared to early-type spirals 
particularly at high redshifts. These are precisely the 
type of systems within which we find fewer bars. More- 
over, in the same figure, the higher mass galaxies also 
have a higher fraction of ordered motions than disor- 
dered motion, although the trend is hard to see in the 
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relatively modest sample size in the high redshift bins. 
These data suggest that the lack of bars may therefore 
be related to the dynamic hotness and the mass surface 
density of these disks. We are currently identifying bars 
and measuring the bar fraction in this sample of galaxies 
and should be able to provide a direct answer for the said 
hypothesis (Sheth et al., in preparation). 

Simulations show that interactions speed up 
bar formation in direct encounte rs, but have lit- 
tle ef fect in retrog r ade ones dToomre fc Toomre 
19721 iNoeuchil [19871: iGerin Combes fc Athanassoulal 



1990; ISteinmetz fc Navarrol 12002ft. in _good agree- 
ment with observations dKormendv &: Norman! 119791 : 
lElmegreen fc Elmegreenl Il982f) . Thus one might have 
expected higher rates of bar formation at high z, where 
interactions are common. On the contrary, it is possible 
for mergings to destroy or seve rely weaken the bar , 
without destroying the disk (e.g. iBerentzen et alj [20031 . 
and references therein). More modeling needs to be done 
before we can say with any certainty what the combined 
effect of interactions and mergings is. Note that we 
discarded from our statistics obviously interacting 
systems based on tidal features or obvious distortions. 
However, if a galaxy is weakly interacting, it would be 
difficult to distinguish it from a non-interacting system; 
this is already the case even in the local Universe. So 
our sample of galaxies is most likely probing quiescent, 
post-merger or weakly interacting disks. 

5.3. The Downsizing Signature in Formation of 
Galactic Structure 



Galaxy "downsizing" was coined bv lCowie et~all(p9fih 
to refer to an evolutionary history in which the most mas- 
sive galaxies formed first. There is strong observational 
evidence for the downsizing phenome na: the presence of 
mass i ve systems at high redshifts (e.g.. lDaddi et al.ll2004L 
I2005t iScarlata et al.ll2006ll2007ft . an orde r of magnitude 
declin e in the typical star formation rate (|Arnouts et al.l 
2005), a change in the star formation activity to lower 
mass systems with decreasing redshifts (|Treu et al.ll2005l; 
Glazebrook et al.ll20"ol iFontana et al.ll2004t iMaier et all 
2005Ll2006HBundv et al.ll2006ft . andthe decrease in char- 



acteris tic luminosity of active gal actic nuclei at low red- 
shifts (|Peilll995tlUeda et al.ll2003ft . 

The results presented here show a downsizing signa- 
ture in the formation of bars. The most massive, lumi- 
nous and red galaxies have a higher bar fraction in the 
highest redshift bin with the most massive and luminous 
systems having a bar fraction close to the present-day 
value. The early presence of bars in these galaxies in 
the context of bar formation ( ^5.21) suggests that these 
systems "matured" early, i.e. they became dynamically 
cool and sufficiently massive to host bars at z> 0.8. In 
contrast, the lower mass systems which are also bluer, 
acquired a majority of their bars at z~0.8. Thus the 
downsizing phenomenon must be considered to be an ef- 
fect more fundamental than one concerned solely with 
the regulation of ongoing star formation; it seems to be 
intimately connected with the dynamical maturity of the 
stellar disk. 

5.4. The Co-Evolution of Bulges & Bars 

Figure [5] shows that in our highest redshift bin there 
is a somewhat higher fraction of bars in galaxies with 



more massive bulges. This is consistent with the struc- 
tural downsizing discussed in the previous section be- 
cause galaxies with bulges are denser and more evolved 
in the center than galaxies without bulges. Bars and 
bulges apparently form at about the same time, with 
later times for lower mass galaxies. Some care is nec- 
essary in the interpretation of this result, however, 
since bulges are an inhomogeneous class of objects. 
In this paper, we defined the bulge as the compo- 
nent in the central part that contributes extra light 
above an extrapolated exponential fit to the outer part. 
This definition includes both classical (3D) bulges and 
disk-like pseudo-bu lges ([Kormendy fc Kennicuttl 120041 : 
lAthanassoulal 120051 ). Pseudo-bulges form by gas in- 
flow and star formation. Because bar s drive inflow 
(jSheth et all 120051 : ISakamoto et~aT] 1 1999ft . there should 
be a correlation between disky bulges and bars, indepen- 
dent of dynamical downsizing. Bars may also contribute 
to the growth of nuclear black holes if they correlate with 
bulges, because there is a tight correlation between black 
hole mas s, stellar velocity dispersion, and luminosity 
of bu l ges dKormendv fc Richstone1[l99 5; Mago rrian et al.l 
119981 : iFerrarese k Merrittl l2000ft . It will be interesting 
to test such differences with subsequent analysis of the 
bulges and bars in the COSMOS data. 

6. CONCLUSIONS 

Bars are an important signpost of galaxy evolution 
because once a galaxy disk is sufficiently massive, dy- 
namically cold and rotationally supported it forms a bar. 
Therefore the evolution of the bar fraction over time is 
an important indicator of the evolutionary history of disk 
galaxies and the assembly of the Hubble sequence. Us- 
ing a detailed analysis of 2,157 L* face-on, spiral galaxies 
from 0.0<z<0.84 in the COSMOS 2-square degree sur- 
vey we have investigated the evolution of the bar fraction 
over the last 7 Gyr. We have undertaken an extensive 
and careful analysis of selection effects (k-correction, sur- 
face brightness dimming, inclination, spatial resolution, 
etc.) which is detailed in the Appendix. Our main results 
are as follows: 

• The bar fraction for L* galaxies drops from about 
65% in the local Universe to about 20% at z=0.84. Over 
this redshift range the fraction of strong bars (SB) drops 
from about 30% to under 10%. Thus at a lookback time 
of 7 Gyr, when the Universe was half its present age, fun- 
damental aspects of Hubble's 'tuning fork' classification 
sequence had not yet fallen into place. Only about one 
fifth of all spiral galaxies were "mature" enough (dynam- 
ically cold, massive and rotationally supported) to host 
galactic structures of the type we see today. 

• For the total f& ar (SB+SAB), the change is far less 
dramatic between z=0.3 and z=0.0 indicating slow evo- 
lution in galactic structures in L* galaxies over the last 4 
Gyr. It is likely that there is significant evolution in the 
formation of bars in the sub-L* galaxies over this period. 

• One of the most significant findings in this study is 
the correlation between % ar and the galaxy mass, lumi- 
nosity and color. We find that in the highest redshift bins 
ibar is higher in the more massive, luminous and redder 
systems. In fact, in the most massive systems, ff, ar is 
already as high at z=0.8 as the local value. These sys- 
tems thus had already arrived with their present Hubble 
types at a lookback time of 7 Gyr. In the subsequent 3 
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Gyr, from z=0.84 to z—0.3, the lower mass, bluer systems 
evolved more slowly toward their present Hubble types. 
Thus the signature of downsizing is intimately connected 
with dynamical maturity of disks and is present in the 
formation of galactic structure. 

• Finally, we find a slight preference for barred galax- 
ies to be more bulge-dominated in the high redshift bin. 
This correlation is consistent with the dynamical down- 
sizing found for bars in general if bars and bulges both 
form earlier and more prominently in the most massive 
galaxies. The lack of a stronger correlation may be re- 
lated to the variety of bulges: bars are also likely to be in- 
volved with the inflow that builds pseudo-bulges. Given 
the strong correlation between bulge properties and black 
hole mass seen today, there may be a co-evolution of bars, 
bulges and black holes in some galaxies. The exact de- 
tails of these processes remain to be investigated. 
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APPENDIX 

Appendices 

ANALYSIS OF SELECTION EFFECTS 

Although we have carefully chosen a robust sample of galaxies, used multiple methods for identifying bars and 
analyzed a sample of local SDSS galaxies in the same manner as the COSMOS galaxies our results are in 

contradiction to some previous studies. Therefore we do additional investigation of the remaining possible selection 
effects (cosmological / surface brightness dimming and spatial resolution), which might produce a declining bar fraction. 

K- Corrections and the Effects of Bandshifting 

The ACS data for COSMOS utilizes the broad F814W filter which traces different rest-frame wavelengths at different 
redshifts. As a result it is imperative to understand the effects of k-correction (bandshifting) and correct them as 
necessary. To quantify the effects of k-correction on the identification of bars, we examined a local sample of 139 
galaxies in all five Sloan bands (u,g,r,i and z). We selected nearby (< 100 Mpc), face-on (b/a > 0.58), large (a 
90% radius > 2 kpc) and bright (Mb -19.7, Mb estimated from g — r color s and g-band mag nitudes iBlanton et al.1 
2003) spiral galaxies from the S loan Digital Sky Survey (SDSS) (|York et al.ll2000t iGunn et al.lll998f) Dat a Release 4 
( Adelman-McCarthv et aLl l2006). The data were mosaicked and calibrated using the methods described by I West et al.l 

We used the same bar identification methods (ellipse fitting & visual classification, Ej3j) for the SDSS data as for the 
COSMOS data and measured the bar fraction in each band. The results are shown in Figure \7\ The bar fraction 
is un changed from the z-band to the g-band at ib ar ~0-6- This is consisten t with a number of previous studies 
fe.g.. iMenendez-Delmestre et al.l [20071: lEskridge et al.ll2002l; IWhvte et al.|[2002l ) that have shown that the overall bar 
fraction does not change appreciably between the optical and near-infrared bands. This figure also demonstrates that 
the shifting rest-wavelength of observation in our sample does not bias our measurements of the bar fraction, provided 
we restrict the maximum redshift of our sample appropriately. 

The maximum redshift chosen is important, because Figure [7] shows that the SDSS bar fraction does appear to decline 
markedly in the u-band. At this wavelength, in a majority of cases, the ellipse fitting technique fails completely. This 
is not unexpected, and a component of this decline may find its origin in the relatively poor signal-to-noise of the 
SDSS u-band data. However, we suspect that the bulk of this decline is real. Bars are primarily stellar structures and 
some become significantly fainter and sometimes disappear altogether short wards of the Balmer break. A dramatic 
example of this is shown for the nearby strongly barred spiral NGC 4303 in Figure 1 of She th et al.l (|2003f ) . This is 
further justification for our chosen limiting redshift in this paper, because by restricting our sample to z — 0.835, the 
F814W filter does not probe bluer than rest-frame g-band, as shown in Figure [H 
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Fig. 7. — Bar fraction as a function of the SDSS filters (u,g,r,i and z) for a local sample of 139 SDSS galaxies. The left panel shows 
the results using the visual classification method and the right panel using the ellipse fitting method described in J3] These are the same 
methods used for the analysis of the COSMOS data. The asterisks show the total bar fraction and the diamonds show the strong bar 
fraction, as described in The pair of numbers above each point are the total number of galaxies and bars identified by each method. 
The u band data point is missing in the right panel because the ellipse fitting algorithm fails in a majority of the galaxies in the u band. 
As noted in the text the main point of this exercise is to quantify the effects of k-correction on bar identification. We find that there is a 
significant k-correction for the bar fraction short wards of the Balmer break in the u-band but the bar fraction is constant from the z-band 
to the g-band. 
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Fig. 8. — The rest-frame wavelengths traced by the F814W filter over the redshift range of this study. At the highest redshift z=0.835, 
the F814W filter is traci ng th e rest-frame SDSS g band filter where the effects of k-correction and identification of bars is still minimally 
affected as discussed in ilA.ll 
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Fig. 9. — As for Figure [T] except that galaxies with clump cluster, compact or other non-standard morphologies are discarded (§3). The 
overall trends noted in the discussion of Figure [l] remain the same. 

Objects with Peculiar Morphology 

We considered the possibility that the bar fraction may be incorrectly measured at high redshifts due to the pres - 
ence of a more exotic variety of morphologies that hav e been observed at high redshifts. lElmegreen et al.l |2007); 
lElmegreen fc Elmegreenl ((20061 . [2005aD ; [Elmcgrce n et al.l (|2005bl lcT) have identified clump-cluster, compact, and chain 
galaxies at high redshifts which do not have counterparts in the local Universe. Could the f& ar . be artificially lower 
because we are including more of these objects in the total sample at the higher redshifts? As noted earlier, this is 
not the case for the visual classification method for which these objects were already excluded. We specifically exclude 
these objects for the ElPa classification method and recompute the evolution of the bar fraction. The results are shown 
in Figure [9l The drop in the bar fraction in this revised sample is less steep as expected. ffc ar drops to 0.28±0.2 at 
z=0.835, and £sb drops to 0.12±0.02. Both methods of classifying bars thus show declines of 50% in the bar fraction 
from the present day to z=0.835. 



Cosmological Surface Brightness Dimming 

Surface brightness dimming is critical even at z< 1 because it evolves so strongly with redshift (oc (1 + z) 4 ). The 
measurement of the bar fraction may be affected if the data are too shallow because as the disk of the galaxy fades 
the bar, which usually has a higher mean surface brightness, may remain visible and therefore be misclassified as an 
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Detectability of Outer Disk Isophotes by different surveys 
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Fig. 10. — A plot of the noise to signal ratio for 0.6"diameter galaxy with a rest-frame Hi a b = 24.5 magnitudes arcsec 2 in COSMOS. 
The choice of 0.6" diameter was based on previous studies of galaxy sizes lIFerguson e t al. 2004). Our own data are consistent with these 
size es timates. We find that the median exponential scale length of L* galaxies in COSMOS is 3.1 kpc (OV 39) at z=0.835 ( Sh cth et al.l 
2007b). In the sample analyzed here, 93% of the galaxies have a half light radius greater than than 0.3". So the noise to signal measurement 
shown here for COSMOS data for the outermost disk isophote is a lower limit. The horizontal lines show the la sensitivity limits for the 
COSMOS, GEMS, GOODS, and UDF survey. Note that GEMS will be adversely affected in measurements f bar even for L* galaxies at 
z > 0.5 because of the low S/N in the underlying disk. 



inclined spiral (e.g. JJogee et alJ l2002). To investigate this possibility we investigated the ability of the COSMOS data 
(and a few other surveys) to trace an outer disk (several times the typical bar semi-major axis) isophote as a function 
of redshift . 

The median bar semi-major axis meas ured from the 2MASS Large Gala xy Atlas survey is a& ar =4.2±2.9 kpc and 

relative S1ZG IS &>l)ar 

I R 25 = 0.29±0.17 ijMenendez-Delmestre et alj|2007l) . where R25 is the classic 25 magnitude 
arcsec -2 isophote in the B-band. Thus the R25 radius is at least three times larger than a typical bar. In the outer 
regions of a galaxy, B-I = 1.5. So the corresponding limiting I-band isophote is at [ii = 23.5. An isophote a full 
magnitude fainter (/i/=24.5 magnitudes arcsec -2 ) can therefore be safely considered to be an outer disk isophote. 

In Figure [POl we show the noise to signal value reached for a /i/=24.5 magnitudes arcsec -2 isophote fading due to 
cosmological surface brightness dimming. The differently colored lines are for the different SED types. It is clear from 
this figure that COSMOS, GOODS, HDF and UDF are sufficiently deep to allow one to detect the outer edges of a 
typical L* disk to z~l. This is not the case for the shallower GEMS data which are unable to quantify the outer disk 
isophotes for L* galaxies at z > 0.5. 

In this calculation we assume that the disks have evolved passively from z=l to the present, when in fact there is 
a magnitude of luminosity evolution in disks from the increased star formation rate. Therefore choosing a luminosity 
evolution in the sample selection criteria, as we have done for the COSMOS data, further improves the signal to noise 
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Fig. 11. — A plot of the total bar fraction (ff, ar ) versus measured surface brightness in the COSMOS sample in three redshift bins. We 
find that there is no correlation between the bar fraction and the surface brightness of the galaxies, as might have been expected if surface 
brightness dimming was affecting the measurement of the bar fraction. The overall trend of lower bar fraction with redshift can be seen 
readily in these panels. The results are the same when considering only the strong bar fraction. 

and ability to confidently trace the outermost disk isophotes with redshift. Based on this analysis we are confident 
that disk dimming is not responsible for the observed decline in the bar fraction with redshift in the COSMOS data. 

We further tested the effects of surface brightness dimming on the data with a second empirical check by comparing 
the bar fraction as a function of the observed surface brightness of the disks in each redshift bin. This test is shown 
in Figure 1111 If indeed surface brightness dimming was responsible for a decline in the bar fraction and we were 
preferentially classifying lower surface brightness disks as unbarred, we should see a correlation of ibar with /i. We 
see no trends in the bar fraction with [i in any redshift bin. Thus we conclude that the observed evolution of the bar 
fraction shown in Figure [T] is not due to cosmological surface brightness dimming of the disks. 



Spatial Resolution and the Bar Fraction 

One possibility for the observed decline in the bar fraction could be th at we are preferen tially missing small bars 
at higher redshifts. As already noted earlier, (and as shown in Figure 3 of iSheth et al.l [20031 ) ACS resolution (1 pixel 
= 0'.'05 = 0.38 kpc at z=0.835) is adequate for identifying all bars larger than 2 kpc at all z<0.835. Bars smaller 
than this, at least in the local Universe, are nuclear bars and not the pri mary bars tha t we are concerned with in 
this paper. Very small bars in late type galaxies have been measured by lErwi n (2005) in a local sample but it is 
unclear whether these galaxies are comparable to the large, L* and brighter galaxies being examined in our sample. 
For L* galaxies, there is substantial evidence that the size of galaxy disks does not change significantly to z^l. For the 
COSMOS sample, the median disk sc ale length is unchanged oyer the redshift range under investigatio n, consistent 
with similar results found previously (jRavindranath et all 12004 iBarden et al.ll200"5tlSargent et al.ll2007h . Given that 



the bar semi-major axis is typically O.3R25, the lack of a change in disk sizes also indicates that it is unlikely that bars 
change their sizes significantly as a function of redshift. 

As another check we decided to see if the measurement of the bar fraction was affected by the size of galaxies in a 
given redshift bin. If the bar fraction was indeed declining because of smaller galaxy disks with unresolved bars, one 
would expect to see a decreasing ff, ar with smaller disks and this effect would be most pronounced at higher redshifts 
where the linear resolution of the ACS data is the coarsest. When we plot ibar versus the half light radius of galaxies in 
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Fig. 12. — A plot of the total bar fraction (f(, ar ) versus half-light radius. If the bar fraction was indeed declining because of smaller 
galaxy disks with unresolved bars, one would expect to see a decreasing fj, ar with smaller disks particularly in the highest redshift bin. We 
find no such trend and therefore conclude that the general decline in the bar fraction is not correlated with disk size. 

different redshift bins in Figure [T2l we find no significant correlation between the bar fraction and the size of the disk in 
any of the redshift bins. There is a slight decline in the bar fraction for the smallest galaxies. This reflects the finding 
that bars are less frequent in less massive systems and not a selection effect. The result is the same when considering 
only isB, or using the exponential scale length for an estimate of the galaxy disk size. We therefore conclude that 
there is not a preferential loss of bars in disks of any particular size at any redshift in this study. 

ARTIFICIALLY REDSHIFTING GALAXIES 

The above sections have addressed most of the classical selection effects that plague high redshift studies, namely k- 
correction, resolution, surface brightness dimming, etc. To conclusively test all of these eff ects we artificially reds hifted 
the g-band images of all 139 SDSS galaxies to z=0.84 following the technique outlined in lGiavalisco et aTl (|1996[ K We 
redshifted the galaxies to the F814W filter, rebinned the image, took into account surface brightness dimming and 
matched the noise characteristics of our ACS data. Examples of the artificially redshifted galaxies are shown in Figure 
IT3l We re-classified these images into strongly barred, weakly barred and unbarred spirals as we had done before. In 
127/139 galaxies the classification remained unchanged. Of the remaining twelve galaxies, seven are now classified as 
weakly barred in these images compared to the original SDSS images where they were classified as unbarred. Four 
galaxies classified previously as weakly barred are now classified as unbarred. And one galaxy that was classified 
as strongly barred is now classified as an unbarred spiral. So overall the bar fraction did not change appreciably 
(original SDSS images f& ar =0.59, artificially redshifted SDSS images fi, ar =0.60). These results are remarkable but not 
unexpected and reflect the exquisite sensitivity of the ACS survey, which was designed to be as deep at z=l as the 
SDSS survey is locally. 
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Fig. 13. — Two examples of artificially redshifted SDSS galaxies. A barred spiral is shown in the top row and a spiral galaxy in the 
bottom row. The left column shows the original SDSS g-band image and the right column is the artificially redshifted image at z=0.84 
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TABLE 1 

Galaxy Classification & Bar Fraction with Redshift 



ZL a z v h N bm c SP Mn d WB[j„ c SB Wn f f bar S la [ f SB h la 1 



Visual Classification 



SDbb-i 


0.00 


139 


57 


23 


59 


0.59 


0.04 


0.42 


0.04 


0.14 


0.26 


75 


30 


19 


26 


0.60 


0.06 


0.35 


0.05 


0.26 


0.37 


243 


108 


63 


72 


0.56 


0.03 


0.30 


0.03 


0.37 


0.49 


254 


136 


46 


72 


0.46 


0.03 


0.28 


0.03 


0.49 


0.61 


265 


161 


52 


52 


0.39 


0.03 


0.20 


0.02 


0.61 


0.72 


504 


336 


78 


90 


0.33 


0.02 


0.18 


0.02 


0.72 


0.84 


364 


251 


50 


63 


0.31 


0.02 


0.17 


0.02 






Ellipse/ PA Class! 


ification 












SDDS-i 


0.00 


139 


57 


26 


56 


0.59 


0.04 


0.41 


0.04 


0.14 


0.26 


83 


29 


32 


22 


0.65 


0.05 


0.27 


0.05 


0.26 


0.37 


267 


114 


84 


69 


0.57 


0.03 


0.26 


0.03 


0.37 


0.49 


282 


134 


85 


63 


0.52 


0.03 


0.22 


0.02 


0.49 


0.61 


326 


210 


68 


48 


0.36 


0.03 


0.15 


0.02 


0.61 


0.72 


668 


495 


97 


76 


0.26 


0.02 


0.11 


0.01 


0.72 


0.84 


530 


412 


71 


47 


0.22 


0.02 


0.09 


0.01 



a z/, : Lower limit of redshift bin 

b Z[j: Upper limit of redshift bin 

c N(, in : Total number of spiral galaxies in bin 

d SPf,i n : Number of galaxies classified as unbarred in bin 

° WB(, in : Number of galaxies classified as weak bars (SAB) in bin 

f SB Mn : Number of galaxies classified as strong bars (SB) in bin 

S har = WB Mn + SB Mn / N Mn 

h fss = SB Mn / N bin 

1 Error assuming binomial statistics. l<T=^/~(f(l — f)/N bin ), where f is i bar or f$B 



